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Abstract LS 5039 is one of the few X-ray binaries detected at VHE, and potentially 
contains a young non-accreting pulsar. The outflow of accelerated particles emitting 
synchrotron emission can be directly mapped with high resolution radio observa- 
tions. The morphology of the radio emission strongly depends on the properties of 
the compact object and on the orbital parameters of the binary system. We present 
VLBA observations of LS 5039 covering an orbital cycle, which show morpholog- 
ical and astrometric variability at mas scales. On the other hand, we discuss the 
possible association of LS 5039 with the supernova remnant SNR G016.8— 01.1. 



1 The binary system LS 5039 



LS 5039 is a high mass X-ray binary containing a compact object of unknown mass 
(1.5-1OM ) that orbits an 06.5 V((f)) star every 3.9 days. The system has an eccen 



trie o rbit (e = 0.35), and it is located at 2.5 kpc dCasares et all 120051 : lAragona et al 



2009). LS 5039 is one of the few known binaries that displays high-energy (HE) 
and very high ener gy (VHE) gamma-ray emission (E > 100 GeV) with a clear 



orbital modulation ( Aharoni an et al. . 2006t Abdo et al. . 20091) . Several theoretical 



models have been developed to explain the multiwavelength orbital behavior of 
LS 5039. The HE/VHE emission is basically interpreted as the result of inverse 
Compton upscattering of stellar UV photons by relativistic electrons. The acceler- 
ation of electrons seems to be produced by shocks between the relativistic wind 
of a young non-accreting pulsar and the w ind of the stellar companion (iDubus . 
2006; Sierpowska-Bartosik & Torres, 2008). No radio pulses have been detected at 
1.4 GHz lMorris et alj (12002b . although the strong free-free absorption with the stel- 
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lar wind may prevent any detection of pulsed radio emission. Therefore, the nature 
of the compact object (black hole or pulsar) is unknown. 

In the non-accreting pulsar scenario the shocked material is contained by the stel- 
lar wind behind the pulsar, producing a 'bow' shaped nebula extending away from 
the stellar companion. The high energy emission is produced in the region where the 
wind pressures are balanced, while a tail of accelerated particles forms behind the 
pulsar. The cooling processes of these accelerated particles along the adiabatically 
expanding flow produce the non-thermal broadband emission from radio to X-rays. 
The morphology is similar to the one expected in isolated pulsars moving through 
the ISM but, as a consequence of the orbital motion of the binary system, the tail of 
the flow bends following an elliptical path during the orbital cycle. 



2 Study of the milliarcsecond radio emission of LS 5039 

Free-free absorption on the stellar wind becomes optically thin at distances ~ 1 AU. 
At scales of a few AU the radio structure produced by the outflow becomes relevant. 
LS 5039 is located at 2.5 kpc, and hence the morphology can be studied at scales of 
1-10 mas (2.5-25 AU), which are directly observable with VLBI. In this context, 
the scenario provides two main predictions about the radio behavior at mas scales. 
On one hand, it is expected that the direction of the extended emission changes 
with the pulsar's orbital motion. On th e other hand, th e peak of the radio emission 
should trace an ellipse of a few mas dDubusl 12006)) . The orbital changes of the 



emission at different frequencies contain relevant information on the nature of the 
compact object, the orbital parameters of the system, the geometry of the region 
where particles are accelerated, and the processes involved. 



3 Radio observations 



The first VLBA observation of LS 5039 at 5 GHz, conducted in 1999, showed a 
^20 mJy central core and bipol ar extended emission of a few mJy extending over 
6 mas on the plane of the sky (Parede set al. , 2000). Similar observations at two 



different orbital phases ( Ribo et al. , 2008b showed a changing morphology at mas 



scales: the core component had a constant flux density, and the Position Angle (PA.) 
of the direction of the elongated emission changed by 12 ± 3° between both runs, 
with a remarkable symmetry change. Therefore, morphological changes at 5 GHz 
occur on timescales of the order of the orbital period. 

To further study the morphological variability, we observed the source in 2007 
during five consecutive days, to cover a full orbit of 3.9 days, and an extra day to 
disentangle between orbital or secular variability. The observations were centered at 
the orbital phases 0.98, 0.23, 0.49, 0.75, and 00, each one spanning 6 hours (0.07 
in phase), computed using the ephemerides of lAragona et alJ (120091) . The observa- 
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tions were phase-referenced to the phase calibrator J1825 — 1718, and an astromet- 
ric check source was observed regularly during the runs. The self-calibrated images 
show a main core and extended emission up to ~ 6 mas, as in previous observations. 
The P.A. of the extended emission from the main core changes significantly every 
day, covering angles between —50 and —85° for phases in the range 0.5-0.0, and 
P.A. 120° at phase 0.23. Therefore, a subtle change in morphology happens in less 
than 20 hours after periastron. 

The resolution of the phase-referenced images is limited by the scatter broaden- 
ing of the calibrator. The astrometric accuracy, extrapolated from the mean disper- 
sion of the astrometric check source, is 0.23 and 0.27 mas in right ascension and 
declination, respectively (Molddn et al., in preparation). The positions of the peak 
of the emission for phases 0.5-0.0 are compatible within errors. The peak of the 
emission at phase 0.23 shows a displacement of 2.3 ±0.3 mas in P.A. 26°, nearly 
opposite to the extended emission sense. This indicates that the variations of the 
extended emission with respect to the core are a combination of intrinsic variability 
of this extended emission and an absolute displacement of the core component. 



4 Proper motion 



As a by-product, these phase-ref erenced obs ervations provide one average precise 
position of the source in the sky. iRibo et al.l (120021) computed the trajectory of LS 
5039 for the last 10 5 yr using optical/radio astrometry from 1905 to 2002. Their 
result marginally suggests an association with SNR G016.8— 01.1. This new VLBA 
position, combined with previous radio interferometric measurements spanning 9 
years, allow us to calculate a more precise proper motion and past trajectory of the 
source (see solid line in Fig. [TJl. The new past trajectory of LS 5039 is compatible 
with the center of the SNR. With a firm association, a kinematical age of the compact 
object would be obtained, which would have direct implications on its properties. 
In case the system contains a pulsar, it would probably be a young pulsar in the 
non-accreting phase, and therefore, compatible with the shocked winds scenario. 



5 Discussion 

The emission at 5 GHz traces the short-lived electrons of the outflow up to 15 AU. 
The morphological changes along the orbit allow us to model the velocity, the en- 
ergy, and the cooling time-scale of the flow of particles that originates the extended 
emission. The morphological and astrometric information constrains the inclination 
of the orbit (/), a key parameter of the system. The mass function of the system and 
i yield the mass of the compact object. On the other hand, it should be possible to 
clearly trace the peak position along the orbit with observations at higher frequen- 
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Fig. 1 Wide field radio 
map of LS 5039, showing 
the position of the nearby 
SNR G016.8-01.1. The grey 
scale emission is at the 6 cm 
wavelength. The overlaid 
contours correspond to the 
NVSS map at the 20 cm. 
The arrows mark the proper 
motion senses. The dashed 
line is the computed trajectory 
reported in Ribo et al. (2002), 
whereas the solid line is the 
computed trajectory from our 
fit (see text). The error boxes 
are calculated for 10 s yr ago. 



cies (where the phase calibrator is much more compact), yielding direct information 
on the absorption around the system. 

The VLBA images from 2007 at 5 GHz contain relevant information to test the 
models (Moldon et al. in prep.), but the lack of continous astrometric information 
following the peak of the emission at all orbital phases is a key point to better con- 
strain the physical properties of the system. If the peak position is expected to be 
shifted between 1-2 mas, it is not possible to unambiguously distinguish the dis- 
placement of the peak from intrinsic variations of the extended emission. The flow 
velocity and cooling times strongly depend on this ambiguity, which can be disen- 
tangled with accurate astrometry, to be obtained in the future. 
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